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We characterize the degree of entanglement of a subsystem of k particles in a A'^-two level system 
(fc < N/2) initially prepared in a mesoscopic superposition |V) = / dO f {e){\<f)i{e))'^^ + \<))2{0)Y'^'^), 
where f{6) is a gaussian or a delta function, subject to the time evolution described by a dephasing 
channel. Negativity is used as a measure of entanglement for such system. For an arbitrary number 
of particles A'^, numerical results are given for the full time evolution up to ten particles. Analytical 
results are obtained for short times and asymptotic time regimes. We show that negativity is 
initially proportional to the square root of the product of the number of particles in each partition, 
the overlap \{<f)\{6)\(f)2(0))'^ and the coupling to the environment. Asymptotically, negativity tends 
to zero, a necessary condition for separability. 
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Generation of macroscopic superpositions has always 
been a subject of interest, since they may be viewed as 
a key to the understanding of some of the striking dif- 
ferences between the quantum and classical descriptions 
of the world such as the "nonexistence at the classical 
level of the majority of the states allowed by Quantum 
Mechanics 0, Q. Superposition of few photons in a 
cavity has been demonstrated in Ref.Q- Recently, the 
observation of macroscopic superpositions has been pro- 
posed for different many particle physical systems, such 
as Bose-Einstein Condensate (BEC)jj|, superconducting 
devices (SQUIDs) "sl and atomic ensembles [6|. Quan- 
tum Interference in a multi-particle system has been in 
fact observed in superconducting devices (SQUIDs) 0. 
All these phenomena are of relevance for Quantum Com- 
putation, especially because quantum computation re- 
lies on the creation and maintenance of complex entan- 
gled states of many subsystems (superpositions of many 
particle systems). However, this is a difficult task be- 
cause of decoherence. The interaction between the sys- 
tem and the environment creates entanglement between 
them and at the same time reduces the entanglement 
within the system itself. An important question in this 
context is how to quantify and classify entanglement. For 
example, a way to compare the effective size of a multi- 
particle macroscopic superposition with the ideal GHZ- 
state (Greenberger-Horne-Zeilinger) has been shown in 
Ref. 0. 

Pure state entanglement of a bipartite system is 
presently well understood On the other hand, entan- 
glement of mixed states is still an open problem in spite 
of all efforts which have been devoted to this problem . 
Things get even more complicated when one is interested 
in investigating the entanglement of multipartite pure or 
mixed states. However, an important step has been given 
in order to determine a computable measure of the en- 
tanglement of any mixed state of an arbitrary bipartite 
system {A^B} This measure is based on the trace 
norm of the partial transpose p^^ of a bipartite mixed 
state p. It measures the degree to which p'^^ fails to be 
positive, and therefore it can be viewed as a quantitative 



version of Peres criterionjllj. This measure is called neg- 
ativity and it is defined as M{p) = X^i I I ; where Vi is 
the ith- negative eigenvalue of the partial transpose p"^-* . 
In addition, negativity is an upper bound to the distil- 
lation of entanglement (an asymptotic distillation rate) 
[To|. one of the fundamental measures of entanglement. 

The aim of the present contribution is to characterize, 
as a function of time, the degree of entanglement of a 
subsystem of k particles in an TV particle system (fc < 
N/2), initially prepared in a mesoscopic superposition 

IV) = I d0f{e){\Me)r'' + 102(0))^^), (1) 

where \M^)) = cos(6')|0) + sin(6')|l) and IMO)) = 
cos(6l)|0) -sin(6i)|l) and f{e) = (27rs2)-i/2 cxp[-(6l - 
Oq)'^ / (2s^)] with s — or s 0, subject to the time evolu- 
tion described by the dephasing channel. The dephasing 
channel is chosen to model decoherence without dissipa- 
tion. The superoperator which describes the evolution 
has the form has the form £{p) = pop -|- (1 — Po)a^.pa^, 
where po = {1 + e~'''*)/2 and az is a Pauli matrix |7lll3|. 
This quantum channel describes the probability of the 
qubit to remain intact or to suffer an " phase flip error" . 
The choice of this particular evolution is related to the 
fact that only in this case an analytical solution for the 
negativity can be found. We consider the case where 
each of the particles is coupled to an independent envi- 
ronment. In that case, the time evolution of a system of 
N particles of two-level swtems (qubits) is given by (see 
also for example 01 Ed) 13) 

p{t) = £{\Mdo)){MOo)\r'' + £{\Mdo)){Mdo)\r'' 
+ £{\MOo)){Meo)\r'' + £i\M0o)){M0o)\r''- 

Negativity: We calculate the time scale of entangle- 
ment of a bipartite {k, N — k} system of N-particles de- 
scribed by Eq. (QJ under the time evolution of the dephas- 
ing channel. As a measure of entanglement, we use the 
negativity A/'(p). In this case, the negative eigenvalues 



2 



of the partial transpose of the density matrix p with re- 
spect to k subsystem have to be evaluated. If the state is 
separable, the negativity A/'(jo) = 0. In ref.0, the nega- 
tivity is used to determine the decay of an upper bound 
of the entanglement of distillation of two optical qubits 
interacting in a lossy nonlinear cavity. To our knowl- 
edge, aside from important mathematical results, there 
are no explicity numerical or analytical results regarding 
the degree of entanglement of a particle state subject 
to a particular time evolution in the literature in which 
negativity is used as a quantitative measure. 

(5-function form for f{d): Let us first consider the 
simplest case 5 = 0, where f(9) becomes a delta function. 
For t = 0, there are — 4 degenerate null eigenvalues 
and four non-vanishing eigenvalues of the partial trans- 
pose p^'' , where one of them is negative, showing that 
the system is entangled at t = 0. These non- vanishing 
eigenvalues {Ai_2,3,4} are given by 
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[2(l + cos(20o)^)] 
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V(l - cos2'=(20o))v/(l - cos2(^-'=)(2eo)) 



[2(l + cos(20o)^)] 
and the corresponding eigenvectors {|Ai_2,3,4)} are 
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with 



- (l'/'i(9o))'^'''±l02(9o))»<'') 



, -, where i represents 

^2[l±cos'(2eo)] 

the number of particles in each partition, i.e., i = fc or 
i^N-k. 

The time evolution raises the degeneracy of the iV — 4 
null eigenvalues, creating N — k new negative eigenval- 
ues if iV — fc > 3. The negative eigenvalues of the 
partial transpose p'^'' (t) can be found for short times 
(7t < 1), by solving the equation p^'=(t)|r) = r(t)|r), 
where r(t)'s are the eigenvalues and |r) — Ci\Xi) + 

^^^^ ^j^;^ d.^'^ Vi'^^ |Ai) are the corresponding eigenvec- 
tors. Some of these new negative eigenvalues are degen- 
erate. For a general partition {fc, N — k}, there are five 
negative eigenvalues where two of them are degenerate. 
From these two negative eigenvalues, one is iV — A: — 1 
degenerate and the other one is fc — 1 degenerate. The 
number of negative eigenvalues is related to the size of 
the system we consider. We see thus that in spite of 
the simplicity of the time evolution (linear and indepen- 
dent coupling to the environment), the appearance of new 
negative eigenvalues is a clear signature of the considered 
many particle system . 

In order to illustrate this effect, we consider the par- 
tition {1, TV — 1} for which analytical results can be ob- 
tained for short times {jt ^ l).In this case, negativity is 
given by 



U{p) = - [(A^ - 2)r3 (t) + Fa {t) + Fi {t)] 



(2) 



where Ti{t), T2{t) and F3(i) are the negative eigenvalues 
and are given by 



r2(t) 
r3(t) 



|A4v/(6+-l)(fe--l) 



with 



A 



B = 



,7^ 



(1 - (iV + 1)^) + 2{N - + (TV + l)^]a_a+ 

+ (|)'[(TV-2)6+ + l][(iV-2)6_ + l] 
4(|)'(1 - (iV + 1)|)'[(A^ - 2)6+ - (iV - l)al + 1] 



[(TV - 2)6+ - (iV - 1)4 + 1] 



(3) 



where 



a± 
b± 



cos(26io)±cos(^-2)(26/p) 

[l±cos(^-i)(20o)] 
cos2(2go)±cos(^-3)(2go) 
[l±cos(^"i)(26io)] 



(4) 



In Fig. (1) we show the time evolution of the negativ- 
ity for short time scales: the full curve corresponds to 
the (numerical) exact solution and it coincides with the 
crosses which show the case where all negative eigenval- 
ues are taken into account (given by Eq.Q). The dotted 
curve corresponds to the evolution of Fi(t). This figure 
illustrates the importance of taking these new eigenval- 
ues in order to calculate the negativity. In addition it 
shows the perfect agreement between the numerical ex- 
act solution for short times and the analytical calculation 
of negativity valid for short times, given by Eq.m. For 
N = 2 and = 3, there is only one negative eigenvalue 
which evolves in time and it is the one which already 
exists at t = 0, A4. 

We are interested in the case of large number of par- 
ticles A^ > 1 where the 6*0 < 1. For the case NO^ < 1, 
one sees that the eigenvalues T2{t) and F3(t) do not con- 
tribute for the negativity. The dominant contribution is 
given by Fi(t). For the partition {1, A^ — 1}, the neg- 
ativity for short times {-ft ^ 1) and 6 1, such that 
NOo-ft ^ 1 , is given by 



(5) 
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FIG. 1: Negativity as a function of time for A*' = 10, 6^0 = tt/S 
and 7=1. The dotted curve corresponds to the evolution of 
Ti{t). The crosses show the case where all negative eigenval- 
ues are taken into account (Eq.l|5J) and it coincides with the 
full curve which corresponds to the numerical exact solution. 

If we consider a general partition {k,N — k}, one sees 
that the structure of the problem does not change dras- 
tically. For this general case in the short time regime 
{-yt ^1), there are five negative eigenvalues, where two 
of them are degenerate. From these two degenerated neg- 
ative eigenvalues, one is — fc — 1 degenerate and the 
other one is fc — 1 degenerate. However, in the case of 
small 6o and ^ 1 , such that N9q <C 1 the main contri- 
bution for the negativity is given by ^lit) (where Ti{t) = 
X4{l + {-tt/2){-N + kA':l:^A^_^'^ + {N^k)Af-'''^A^_^-'''^)), 

with = ((cos(2^?o)±cos(20o)^''"'^)/(l±cos(2^o)**^)) 
where i represents the number of particles in each parti- 
tion k OT N — k). In this limit, the negativity is 

J\f{p) - Vky'iN - k)el{l ~ 2jt). 

This shows that the maximum entanglement is obtained 
for k = N/2, increasing thus the characteristic time 
scale of the negativity. In fact even for a more gen- 
eral evolution: depolarizing channel {£{p) = pop + 
Pi{(TxPcrx + cFyPc^v + cTzpa-z), where po = (1 + 3e"''*)/4, 
Pi — (1 — e~'''*)/4 and a(^x.y,z) are Pauli matrices), we can 
see numerically that one also gets the same dependence 
for the negativity on the number of particles. This is due 
to the fact that the characteristic time scale in this case 
is essentially governed by initial state correlations. 

The full time evolution of the entanglement of the sub- 
system under investigation is given by numerical calcu- 
lations up to ten particles for an arbitrary ^o- Fig- 2 
illustrates the behavior of the negativity for = 10 and 
^0 = 7''/3 as a function of time (numerical exact calcula- 
tion). 

In the asymptotic limit, it is also possible to calculate 
the negativity and it is zero. This is a necessary condition 
for the separability of the system under investigation. 

The particular case 6*0 = 7r/4 in Eq.QJ (GHZ-state of 
the type lip) = |0)f ^ + |1)®^) and the evolution under 
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FIG. 2: The full curve represents the exact numerical solution 
for the time evolution of the negativity for all times. For short 
times, the behavior of this curve can be seen in Fig. 1 and it 
coincides with the analytical solution of negativity given by 
Eq.Q. In the asymptotic limit, the negativity goes to zero. 
In this figure, we consider {A'^ — 1, 1} partition, 9o = ir/3 and 
iV = 10. 



the dephasing channel, the negativity becomes indepen- 
dent of N. One recovers the iV-dependence if \0)x — > 10)^ 
and \l)z- These particular situations illustrate 

the dependence of the time evolution on the initial con- 
dition: If the evolution operator is "orthogonal" (in the 
sense just described) to the initial condition, one slows 
down the entanglement process. Also if the initial state 
of Eq.lU has |0i) = |0) and \(j}2) = cos(6io)|O)-|-sin(0o)|l), 
the negativity is completely governed by the time evolu- 
tion of the negative eigenvalue which exists for t = 0. 
In this case, the state |</>i)(|(/)i| is an eigenstate of the 
evolution operator and we observe that the degeneracy 
among the — 4 null initially existing eigenvalues is not 
removed. 

Gaussian form for f{9): Next we investigate the 
general situation where a finite width s in f{9o) appears 
in Eq.QJ. The gaussian form does not affect the struc- 
ture of the negative eigenvalues. For a general partition 
{A:, A^ — k} and short time regime {'jt ^ 1), there are five 
negative eigenvalues, where two of them are degenerate. 
From these two negative eigenvalues one is A^ — A: — 1 
degenerate and the other one k — 1 degenerate. For the 
partition {1, A^ — 1}, the general form of the eigenvalues 
of Eq. (PI remains unchanged, only the coefficients a± and 
b± are more complicated. Analytical expressions can be 
obtained for cr^ ^ 1 and for short times jt <^ 1. Their 
explicit form will not be given here, since they are not 
particularly illuminating. We see from the numerical im- 
plementation of the analytical result that the negativity 
has the same dependence in the number of particles as in 
the previous case [s — 0). As can be seen from Fig. 3 the 
present situation introduces a slight modification of the 
previous result. Fig. 3 shows the negativity as a function 
of the number of particles A^ for a {A^ — 1,1} partition 
when f{9) is a gaussian or delta function. Fig. 3 was 
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generated from the analytical expressions of negativity 
for short times in both cases (delta and gaussian). 
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FIG. 3: Negativity as a function of the number of particles 
iV, for fixed Oo = 0.1, -yt = 0.02, for {N ~ 1,1} partition. 
The dashed curve represents the negativity as a function of 
the number of particles for the delta case (s — 0) while the 
full curve shows the negativity as function of the number of 
particles for the gaussian case with s = 0.05. 

In summary, we have determined the degree of entan- 
glement of a subsystem of k particles in this N particle 
{k < N/2) system. We observed that negativity is pro- 



portional to the square root of the number of particles in 
each partition, the overlap between the states involved in 
the initial mesoscopic superposition and the coupling to 
the environment , in the limit of small angle and short 
times. In addition, the value obtained here for the nega- 
tivity gives an upper bound for the distillation rate, since 
it has been shown in ref. the relation between these 
two quantities. It would be interesting to investigate 
other dynamical situations, in particular, where the evo- 
lution can be derived from first principles. The case of a 
competition between collective and independent particle- 
environment coupling also deserves further investigation 
in order to see how universal are the conclusions we came 
to. 



Acknowledgments 

ANS acknowledges J. I. Cirac for stimulating discus- 
sions and for suggesting the problem, M. Wolf, J.J. Gar- 
cia RipoU, M. C. Nemes, A. F. R. de Toledo Piza, H. 
Weidenmtiller and M. Weidemiiller for valuable discus- 
sions. 



[1] Letter from Albert Einstein To Max Born in 1954, cited 
by E. Joos in New Techmqes and Ideas in Quantum Mea- 
surement Theory, edited by D. M. Greenberger (New 
York Academy of Science, New York, 1986). 

[2] E. Schrodinger, Die Naturwissenschaften 23, 807 (1935). 

[3] M. Brune et al, Phys. Rev. Lett. 77, 4887 (1996). 

[4] J. I. Cirac, M. Lewenstein, K. Molmer and P. ZoUer, 
Phys. Rev. A 57, 1208 (1998); J. Ruostekoski, M. J. 
Collect, R. Graham and D. F. Walls, Phys. Rev. A 57, 
511 (1998). 

[5] J. R. Friedman et al, Nature 406, 43 (2000); C. H. van 

der Wal et al, Science 290, 773 (2000). 
[6] S. Massar and E. S. Polzik, Phys. Rev. Lett. 91, 060401 

(2003). 

[7] W. Diir, C. Simon and J. I. Cirac, Phys. Rev. Lett. 89, 
210402 (2002). 

[8] C. H. Bennett, H. J. Bernstein, S. Popescu and B. Schu- 
macher, Phys. Rev. A 53, 2046 (1996). 
[9] C. H. Bennett, D. P. DiVicenzo, J. A. Smolin and W. 



K. Wootters, Phys. Rev. A 54, 3824 (1996); P. M. Hay- 
dern, M. Horodecki and B. Terhal, quant-ph/0008134 
V. Vedral and M. Plenio, Phys. Rev. A 57, 1619 (1998); 
L. Henderson and V. Vedral, Phys. Rev. Lett. 84, 2263 

(2000) ; M. Horodecki, P. Horodecki and R. Horodecki, 
Phys. Rev. Lett. 84, 2014 (2000); W. K. Wootters, Phys. 
Rev. Lett. 80, 2245 (1998); G. Vidal, Phys. Rev. A 62, 
062315 (2000). 

[10] G. Vidal and R. F. Werner, Phys. Rev. A 65, 032314 
(2002). 

[11] A. Peres, Phys. Rev. Lett. 77, 1413 (1996). 

[12] J. Preskill, "Lectures Notes for Physics 229: Quantum 

Information and Quantum Computation" (1998). 
[13] B. Schumacher, Phys. Rev. A 54, 2614 (1996). 
[14] C. Simon and J. Kempe, Phys. Rev. 65, 052327 (2002). 
[15] Adam Miranowicz, quant-ph/0402025j 

(2001) . 



